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Abstract 

Nonequilibrium phase transition properties of the +J Ising model under a time dependent oscillating perturbation are 
investigated within the framework of effective field theory for a two-dimensional square lattice. After a detailed anal- 
ysis, it is found that the studied system exhibits unusual and interesting behaviors such as reentrant phenomena, and 
the competition between ferromagnetic and antiferromagnetic exchange interactions gives rise to destruct the dynamic 
first order phase transitions as well as dynamic tricritical points. Furthermore, according to Neel nomenclature, the 
magnetization profiles have been found to obey Q-type, L-type and P-type classification schemes under certain con- 
ditions. Finally, it is observed that the treatment of critical percolation with applied field amplitude strongly depends 
upon the frequency of time varying external field. 
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1. Introduction 

Investigation of disorder effect problems which may 
be originated from random interactions between spins 
or from a random dilution of the magnetic ions with 
non-magnetic species on the magnetic materials has 
a long history, and there have been a great many 
studies focused on disordered magnetic materials with 
quenched randomness where the random variables of a 
magnetic system may not change its value over time. 
Actually, all magnetic materials have some small de- 
fects, and magnetic properties, i.e phase transition tem- 
perature point that separates the ordered phase from 
disordered phase, of sample varies significantly de- 
pending on the type of defects. From the experi- 
mental point of view, a number of studies including 
the site or bond randomness have been devoted to 
better understanding of magnetic properties of differ- 
ent types of magnetically interacting systems such as 
Rb2MneMgi_eF4 m, MnpZni_pF2 A FepMgi^pClj 
(l, Cdi_pMnpTe [4], MneZm^eFz [5], Co(SpSei_p)2 
0, and CopZni_p(C5H5NO)6(C104)2 IJ. 

From the theoretical point of view, a great deal of 
studies have been performed regarding the magnetic 
properties of disordered magnetic systems, and based 
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on the investigation of static or dynamic phase tran- 
sition (DPT) properties of such disordered magnetic 
system, theoretical works can be classified in two ba- 
sic categories. In the former group, static or equi- 
librium properties of these type of systems have been 
analyzed within the several frameworks such as se- 
ries expansion method jS H, [iS \vh, Monte Carlo 
method 112. 13, 14. 15j |16[]. renormalization group the- 




HlMlzinzl S H, repUca method 



efi^ective field theory |l28l |29 



33MM.37, 38, 39] andbethe method lE 



The statistical mechanics of nonequilibrium systems 
is a less well developed and understood field of study 
than that for equilibrium systems. Even though the 
physical investigations regarding these systems bring 
about a lot of mathematical difficulties, the nonequilib- 
rium systems are in the focus of scientists because they 
have an exotic, unusual and interesting dynamic behav- 
iors. For instance, the universality class of the kinetic 
Ising model under a time dependent oscillating mag- 
netic field is different from its equilibrium counterparts 
1 41, i^l- It is exciting to say that a magnetic system 
composed of interacting magnetic moments under the 
influence of a time dependent magnetic field exhibits 
two important phenomena: DPT and dynamic hystere- 
sis behavior Dynamic evolution of the system strongly 
depends upon the amplitude and frequency of the ap- 
plied external magnetic field as well as other intrinsic or 
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extrinsic Hamiltonian parameters. It is possible to men- 
tion that nonequilibrium phase transitions originate due 
to a competition between time scales of the relaxation 
time of the system and oscillating period of the exter- 
nal field. As far as we know, the ferromagnetic system 
exists in dynamically disordered phase at high tempera- 
tures and for high amplitudes of the periodically varying 
magnetic field. In this region, the time dependent mag- 
netization is able to follow the external magnetic field 
with some delay whereas this is not the case for low 
temperatures and small magnetic field amplitudes. The 
mechanism shortly described above points out the exis- 
tence a DPT 14 in . and theoretical point of view, a num- 
ber of studies concerning the DPTs as well as hystere- 
sis properties of different type systems under the time 
dependent perturbation have been performed by using 
well known methods such as Monte Carlo simulation 
1 51 H J4I S IS H S IsSleffective field theory 
I U SEl, mean field theory^ S 13 Q as well 
as hard spin mean field theory [56]. 

In addition to the theoretical published works men- 
tioned above, the DPTs and also hysteresis behaviors 
can be observed experimentally due to recent develop- 
ments in experimental techniques. It is beneficial to 
give a few examples regarding these properties such 
as epitaxial Fe/GaAs(001) and Fe/InAs(001) ultrathin 
films L57J, finemet thin films L58J, [Co(4A°)/Pt(7A°)] 
multi-layer system with strong perpendicular anisotropy 
1I59I1 , Co films on a Cu (001) surface [60], thin poly- 
crystalline NigoFeao films [61J, Fen42Zno.58F2 [62], epi- 
taxial Fe/GaAs(001) thin films ll63ll . epitaxial single 
ferromagnetic fee NiFe(OOl), fee Co(OOl), a nd fee 
NiFe/Cu/Co(001) layers as well as ultrathin ferro- 
magnetic Fe/Au(001) films ll65ll . Based upon the de- 
tailed experimental investigations, it has been discov- 
ered that experimental nonequilibrium dynamics of con- 
sidered real magnetic systems strongly resemble the dy- 
namic behavior predicted from theoretical calculations 
of a kinetic Ising model. From this point of view, it is 
possible to see that there exists an impressive evidence 
of qualitative consistency between theoretical and ex- 
perimental investigations. 

On the other hand, in the latter group there exists 
a limited number of nonequilibrium studies concern- 
ing randomness eff'ects. For instance, by making use 
of both mean field theory and Monte Carlo simulation 
method, it has been shown that the hysteresis loop area 
of a ferromagnet including impurities under a time de- 
pendent magnetic field is a power law function of the 
linear driving rate as A - Aq oc where, Ao, h and y6 
are the static hysteresis loop area, the linear driving rate 
and scaling exponent of the system, respectively ll66ll . 



Very recently, DPT and also hysteretic properties such 
as remanence, coercivity and loop area of the random 
quenched site kinetic Ising model have been probed 
based on the efifective field theory with single-site cor- 
relations and it is observed that the coexistence region, 
where dynamically ordered and disordered phase coex- 
ist depending on selected Hamiltonian parameters, dis- 
appear for sufficiently weak dilution of lattice sites. It 
is also propounds that the considered disordered sys- 
tem indicates the existence of essentially three, partic- 
ularly four types of dispersion curves ||67[ 16811 . Follow- 
ing the same methodology, quenched random bond di- 
luted Ising model, where the spin-spin exchange inter- 
action has a probability p and I - p of taking on values 
J and 0, respectively, has been analyzed by benefiting 
from Glauber type stochastic model lIHoll . and particu- 
lar attention has been devoted the better understanding 
of efifects of bond randomness on the global phase di- 
agrams constructed in related planes and on the micro- 
scopic origin of the magnetic system ItoI After 



a complete detailed analysis, the studied system shows 
that the first order phase transitions as well as dynamic 
tricritical points (DTCPs) disappear while the reentrant 
phenomena exists dependin g on the value of quenched 
bond dilution parameters p 17011 . It is also emphasized 
that frequency dispersions of hysteresis loop area, rema- 
nence and coercivity strongly depend on the quenched 
bond randomness, as well as applied field amplitude and 
oscillation frequency [71]. Furthermore, it is beneficial 
to emphasize that Monte Carlo simulation studies con- 
cerning the nonequiUbrium behaviors and universality 
aspects of diff'erent type Ising spin glasses systems take 
part in literature [72, 7Ji[74l[75||. Even though there ex- 
ists a limited number of studies including the disorder 
eff'ects under the time dependent oscillating magnetic 
field, the dynamic nature of the +J Ising model which 
would exhibit an exotic and unusual dynamic behaviors 
has not yet been investigated. From this point of view, 
in this work, we intend to probe the effects of the ran- 
dom bond dilution process on the kinetic Ising model 
in the presence of a time-dependent oscillating external 
magnetic field by using the effective field theory with 
correlations based on the exact Van der Waerden iden- 
tity for a spin- 1/2 system. 



The remainder of the study is as follows: The dy- 
namic equation of motion and dynamic order parameter 
(DOP) of the kinetic +J Ising model are described in 
Sect. |2] The numerical results and related discussions 
are given in Sect. [3] and finally Sect. |4] contains our 
conclusions. 
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2. Formulation 

We consider a two dimensional +J Ising model de- 
fined on a square lattice, and the considered system is 
subjected to a periodically oscillating magnetic field. 
The time dependent Hamiltonian describing our model 
of magnetic system can be written as 



(1) 



<'■» 



where St - + 1 is the Ising spin variables and first sum in 
Eq. (1) is over the nearest neighbor pairs of spins. We 
assume that the nearest neighbor interactions are ran- 
domly diluted on the lattice according to the probability 
distribution function 

PiJij) = p6(Jij - 7) + (1 - p)6(Jij + J). (2) 

The time dependent magnetic field is as following 



Hi{t) = h cos((L)t) 



(3) 



here, h and w are amplitude and angular frequency of 
the external field, respectively. In order to describe 
the dynamical evolution of the system, we follow a 
Glauber-type stochastic process f69j at a rate of 1/t 
which represents the transitions per unit time, and the 
dynamical equation of motion can be obtained by using 
the master equation as follows: 



d{Si) 
dt 



tanh 



Ei + Hi(t) 



knT 



(4) 



where Ei - Y^j JijS j is the local field acting on the lat- 
tice /, and ks and T denote the Boltzmann constant and 
temperature, respectively. If we apply the differential 
operator technique 11761 17711 in Eq. (4) by taking into an 
account the random configurational averages we get 



dm 
~dt 



— —m + 



SjBi^ F(x)|,=o (5) 



here, where A,y = cosh(y,yV), Bjj - sinh(y,jV) and 
= {(Si)) indicates the average magnetization. V = ^ 
is the differential operator, and the inner (■ ■ ■) and the 
outer (■ ■ ■)r brackets represent thermal and configura- 
tional averages, respectively. When the right-hand side 
of Eq. (|5]l is expanded, the multispin correlation func- 
tions appear. The simplest approximation, and one of 
the most frequently adopted is to decouple these corre- 
lations according to 



{{SiSj---Sl)), = {{Si))r{{Sj))r---{{Sl))r, 



for i j ■ ■ ■ ^ I iItsIi . By making use of Eq. (|6|l, the 
right hand side of Eq. ^ can be easily expanded and 
then, the following dynamic effective field equation of 
motion for the magnetization of the ±J Ising model can 
be found as 



q=4 

dm sr^ I 

— — —m + > £,m 
dt ^ 

,={) 



(7) 



where the coefficients (j can be readily determined by 
employing the mathematical relation exp(Q'V)F(x) - 
F(x + a). Eq. (|7]i is a kind of initial value problem, 
and generally, the solutions act a stable function after a 
certain transition regime. In addition to these, the stud- 
ied magnetic system exhibit two type solutions, and the 
first one is called symmetric solution that obeys the fol- 
lowing property: 



m(f) = —m{t + jt/cj). 



(8) 



In this type of solution, the time dependent magnetiza- 
tion oscillates around zero value which corresponds to 
dynamically paramagnetic phase (P). The second one 
is asymmetric solution where the magnetization oscil- 
lates around a nonzero value which indicates the exis- 
tence a dynamically ferromagnetic phase (F). Accord- 
ing to our solutions, the observed behavior of the mag- 
netization is independent of the choice of the initial 
value of magnetization. On the other hand, depend- 
ing on Hamiltonian parameters and also initial value of 
time dependent magnetization, there exist coexistence 
regions (Fh-P phases) in the phase diagrams in tempera- 
ture versus field amplitude plane. 

The time average magnetization over a full cycle of 
the external magnetic field acts as DOP which is defined 
as follows ilil 



Q = ^^m(t)dt 
where m(t) periodic and stable function. 



(9) 



(6) 



3. Results and Discussion 

In this section, we will focus our attention on the DPT 
properties of ±J Ising model under a time dependent 
magnetic field source. Here, we discuss how the bond 
randomness affects the dynamic evolution of studied 
system in the vicinity of dynamically order and disorder 
transition temperature. On the other hand, it is known 
fact that it is not possible to obtain the free energy for 
kinetic models in the presence of time-dependent ex- 
ternal fields. Hence, in order to determine the type of 
DPT (first or second order), it is convenient to check the 
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temperature dependence of DOR Namely, if the DOP 
decreases continuously to zero in the vicinity of crit- 
ical temperature, this transition is classified as second 
order, whereas if it vanishes discontinuously, then the 
transition is assumed to be first order Keeping in this 
mind, and after detailed analysis with a great success 
to construct dynamic phase boundary (DPB) that sep- 
arates the dynamic ordered phases from the dynamic 
disordered phases, it is shown that the behavior of the 
system changes dramatically with varying Hamiltonian 
parameters. In order to investigate these effects on the 
dynamic behavior of the considered system, we plot the 
DPB in various planes. 

Fig. [T](a-d) represents the DPT lines in reduced tem- 
perature and applied field amplitude plane (ksTc/J - 
hi J) for varying p values and selected applied field fre- 
quencies such as w = 0.25(fl), 0.5(/7), 1.0(c) and n{d). 
Here, the full and dotted lines indicate the continuous 
and discontinuous phase transition points, respectively, 
and also the solid circle shows the DTCP It is observed 
that a reentrant behavior of second order where the two 
successive second order phase transitions appears in the 
relatively low applied field amplitude and low tempera- 
ture regions for low p values. Another important result 
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Figure 1: Phase diagrams of tlie kinetic ±J Ising model in the 
(kaTc/J-h/J) plane for selected bond randomness values with (a) 
Ll) = 0.25, (b) to = 0.5, (c) w = 1.0 and (d) cj = tt. Solid (dotted) lines 
correspond to second- (first-) order DPTs, and • symbols represent 
DTCPs. 



of our study is that an increase in the value of w leads to 
annihilate the reentrant behavior and also to expand the 



dynamically ferromagnetic region. On the other hand, 
with the increasing value of p the dynamic evolution of 
the +J Ising system begins to resemble the pure kinetic 
Ising model, and a dynamically first order reentrant be- 
havior with a coexistence region (F-i-P) phase appears 
where the critical properties of the system depends on 
the initial value of the magnetization. Namely, if ini- 
tial value of the m{t) is selected to be zero, the time 
dependent magnetization oscillates around zero value 
which corresponds to P whereas if it is considered to be 
nonzero, in this time, m(t) oscillates around a nonzero 
value, and this type of behavior corresponds to F. It is 
possible to make an inference that the increasing field 
frequency causes a growing phase delay between the 
magnetization and magnetic field and this makes the oc- 
currence of the DPT difficult, as a result of this mecha- 
nism the DPB gets wider and also the reentrant behavior 
disappears at the low temperature and low amplitude re- 
gions. 




Figure 2: Phase diagrams of the kinetic ±J Ising model in the 
(p - keTc/J) plane for varying reduced magnetic field values with (a) 
CJ = 0.25, (b) CJ = 0.5, (c) w = 1.0 and (d) cj = n. 



Fig. |2] (a-d) illustrates the DPT lines in p and re- 
duced temperature plane (p - kBTdJ) for varying hjj 
values and selected applied field frequencies such as 
(jj - 0.25(fl), 0.5(^7), 1.0(c) and 7r(t/). One can easily 
see that as p value decreases then the ferromagnetic 
phase region gets narrower According to our model, as 
the concentration of active ferromagnetic bonds is de- 
creased from one then the antiferromagnetic exchange 
interactions appear, and due to the competition between 



ferromagnetic and antiferromagnetic exchange interac- 
tions, energy contribution which comes from spin-spin 
interactions gets smaller Hence, the system can un- 
dergo a DPT at lower critical temperatures, and also the 
dynamically ferromagnetic regions get narrower On the 
other hand, the case of h/J - 0.0 corresponds to static 
+J Ising model, in which there is no energy contribu- 
tion originating from magnetic field. Our investigations 
calculated for this applied field amplitude indicates that 
the critical percolation value is pc = 5/6. At the critical 
percolation value the concentration of antiferromagnetic 
bonds equal to 1/6 = 0.167, and it is worth noting that 
this value is in good agree ment with previously pub- 
Ushed works ||29nHlH ]l4l[3ll40ll . Another important 
finding is that increasing h/J values reduce the dynamic 
critical temperature whereas as oj increases then the F 
region gets wider in the related plane. 
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Figure 3: Thermal variations of the DOP curves as functions of the 
selected Hamiltonian parameters h/J, at and p where the letters on the 
curves denote the value of the corresponding Hamiltonian parameter. 



In Figs. Oa-d), in order to elucidate the effects of 
the quenched random bond process on the studied sys- 
tem in detail, we give the thermal variations of DOP 
corresponding to DPBs depicted in Fig. [T]and Fig. |2] 
for a combination of Hamiltonian parameters. Based on 
the upper left and right panels in Fig. |3] it can be eas- 
ily said that as the active ferromagnetic bond concen- 
tration decreases then the system tends to exhibit reen- 
trant behavior of second order Here, two successive 
DPTs exist, and the first one occurs from a disordered 
phase to an ordered phase at relatively low temperature 



regions, while the other one takes place from an ordered 
phase to a disordered phase at higher temperature re- 
gions for decreasing p values. It is also possible to men- 
tion that the aforementioned situation is valid for both 
weak and strong frequency (or amplitude) values which 
can be clearly seen from the upper left and right panels 
in Fig. [3] One of the most interesting findings is that the 
competition between ferromagnetic and antiferromag- 
netic bonds gives rise to the existence of different type 
magnetization profiles. According to Neel nomencla- 
ture 1 8^ 81], it is possible to classify the thermal vari- 
ation of the magnetization curves in certain categories. 
Based on this classification scheme, one can see from 
Fig. [3lb) that considered system exhibits Q-type, L- 
type and P-type magnetization behaviors depending on 
Hamiltonian parameters. Furthermore, the variations of 
the DOP with temperature curves as a function of the 
applied field frequency are shown in the lower left panel 
in Fig. [3]for p = 0.86 and h/J - 0.25. As shown in this 
figure, value of the DOP decreases when the frequency 
approaches the static case. The physical backgrounds 
underlying the behaviors found in Fig. |3lc) are identi- 
cal to those emphasized in above discussions. There- 
fore, we will not discuss these interpretations here. In 
addition to these, the thermal variations of the DOP at 
various h/J values show that depending on the applied 
field amplitudes the system exhibits second order reen- 
trant phenomena in Fig. [3ld). 
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Figure 4: Variation of the bond percolation threshold pc as a function 
of the external magnetic field amplitude h/J for some selected values 
of the oscillation frequency w. The letters on each curve denote the 
frequency value of the external field 



In the following analysis, let us investigate the varia- 
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tion of the bond percolation threshold pc as a function 
of applied field amplitude h/J with some selected val- 
ues of w which is depicted in Fig. |4] We found that 
at relatively small oscillation frequency values such as 
oj - 0.25,0.5 and 1.0, value of pc increases gradually 
then exhibits a plateau which gets wider as oj increases. 
Another important findings related to the percolation in- 
vestigation is that after a certain h/J(- 2.0) value, pc 
value increases rapidly and saturates at pc = 1.0. At 
this point, we should indicate that for p < pc, there 
is no infinite cluster causing to ordered phase and no 
phase transition whereas there is an infinite cluster for 
p > Pc and critical temperature rises from zero. On the 
other hand, it is possible to say that the critical perco- 
lation value Pc strongly depends on a kind of compe- 
tition effect which originates from the collaboration of 
the ferromagnetic and antiferromagnetic bond interac- 
tions with the applied field frequency against the ampli- 
tude of the external field. Moreover, as we mentioned 
above, the case of h/J = 0.0 corresponds to static +J 
Ising model, and the numerical calculations show that 
percolation value is pc - 5/6(= 0.833) which is inde- 
pendent of CO. 
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Figure 5: Bond randomness dependencies of the DOP for fixed tem- 
perature and applied field amplitude at various frequencies of field 
(a). Magnetic field variation of the DOP for considered Hamiltonian 
parameters (b). 



The bond randomness and external field amplitude 
dependence of the DOP are shown in Fig. |5] (a-b) 
for a considered combination of Hamiltonian parame- 
ters. It is obvious from the Fig. |5](a) that an incre- 
ment on the active antiferromagnetic bond concentra- 
tions {I - p) tends to destruct the dynamically ordered 



phase for value of keT/J = 10"^ and h/J = 0.25. We 
notice that the aforementioned situation is dependent on 
applied field frequency, and with the decreasing value of 
CO, the dynamically order-disorder transition point shifts 
to lower value of active antiferromagnetic bond concen- 
trations. At this point, we can also mention that there 
exists a good agreement in a formal manner between our 
Fig. ISfa) and Fig. 5 for a - -I of Ref. (3E] where static 
properties of a quenched random bond Ising ferromag- 
net with anisotropic coupling constants are discussed. 
Moreover, in Fig. |5](b), it is observed that the compe- 
tition between the ferromagnetic and antiferromagnetic 
exchange interactions give rise to destruct the dynami- 
cally first order phase transition as well as to disappear 
the DTCP at various p values and for fixed reduced tem- 
perature and field frequency such as ksT/J - 0.5 and 
0.5, respectively. 
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Figure 6: Variations of DTCP coordinates (a) h /J and (b) 
as a function of bond randomness (p) for considered val- 
ues of oscillating field frequencies such as ai = 0.5 and 1.0. 



As a final investigation, in Figs. |6fa) and|6jb), varia- 
tions of DTCP coordinates h°TCP/j and IcbTDTcp/j with 
respect to the bond randomness are plotted for selected 
oscillating field frequencies such as w = 0.5 and 1.0, 
respectively. As seen in these figures, both the h'-'^'--^/J 
and kfiT'^^^^/J diminish with decreasing active ferro- 
magnetic bond concentration, and the DTCP disappears 
after a certain value. It is also possible to notice that the 
coexistence region in the phase diagrams gets narrower 
and disappear with decreasing p values. Besides, the 
coordinates of DTCPs explicitly depend on the applied 
field frequency, and one can easily deduce from the fig- 
ures that the ordered phase regions and also positions 
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of the DTCPs in related planes tend to expand because 
increasing field frequency gives rise to a growing phase 
delay between the magnetization and magnetic field and 
this makes the occurrence of the DPT difficult. 

4. Concluding Remarks 

In conclusion, we have looked for an answer how 
the time dependent oscillating magnetic field affects the 
DPT properties of ±J Ising model within the framework 
of effective field theory on a two dimensional square lat- 
tice. In our studied model the spin-spin exchange inter- 
action has a probability p and 1 - /? of taking on values 
-1-7 and - J, respectively. For this purpose, the Glauber 
type stochastic process has been used with a great suc- 
cess to determine the time evolution of competing the 
magnetic system. After a detailed analytical and nu- 
merical operations, the DPB separating the dynamically 
ordered and disordered phases has been obtained in dif- 
ferent planes by benefiting from thermal variations of 
DOPs at various Hamiltonian parameters. It is found 
that the competition between +J interactions causes the 
reentrant phenomena for some certain values of amph- 
tude and frequency of the external applied field, and the 
ferromagnetic phase regions get expanded with decreas- 
ing amplitude which is more evident at low frequencies. 
One of the most important findings is that the bond ran- 
domness leads to destruct the dynamically first order 
phase transition as well as to disappear the coexistence 
regions. Namely, after a certain value of ferromagnetic 
or antiferromagnetic concentration, the first order phase 
transitions turn into the second order phase transitions 
and consequently, DTCPs disappear for all frequency 
values. Furthermore, it is observed according to Neel 
nomenclature |80,i81|j that, the magnetization curves of 
kinetic +J Ising model have been found to obey Q-type, 
L-type and P-type classification schemes under certain 
conditions (see Fig. [^b)). We should also emphasize 
that the DPBs in related planes have been constructed 
by making use of the temperature variations of DOP, 
and our calculated phase boundaries do not include spin 
glass phase because there is no infinite cluster causing 
to ordered phase for p < pc, in this region the DOPs 
only indicate the existence of dynamically disordered 
phases. We think that focusing on hysteresis and loop 
area properties of the studied system may be good and 
powerful ways to determine and investigate of the spin 
glass behavior Recently, this type of calculation has 
been done and spin glass phase properties of the d - 3 
+J Ising model have been elucidated in detail within the 
framework of hard spin mean field theory fl56fl . 



On the other hand, it is a well known fact that ef- 
fective field theory takes the standard mean field pre- 
dictions one step forward by taking into account the 
single-site correlations which means that the thermal 
fluctuations are partially considered within the frame- 
work of effective field theory. It is possible to mention 
that the effective field theory can be successfully ap- 
plied to such non-equilibrium systems in the presence of 
bond randomness, however, the true nature of the phys- 
ical facts underlying the origin of the coexistence phase 
and also dynamic first order transitions may be further 
understood by benefiting from a well defined powerful 
method such as Monte Carlo simulation technique. 
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